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Abstract. For the discrete Schrodinger operator we obtain sharp 
estimates for the number of negative eigenvalues. 



1. Introduction 

We study the estimates of the number of negative eigenvalues of the 
discrete Schrodinger operator 

(1.1) H^v = -^-aV 

in the Hilbert space l2{^^\ d > 3. We use the standard notation 
for the lattice points: x = {xi, . . . ,Xd), with Xj e Z for each j. The 
discrete Laplacian is 

d 

(Am)(x) = ^(u(x + Ij) + u{x - Ij) - 2u{x)), X e Z'^, 
i=i 

where Ij is the multi-index with all zero entries except 1 in the position 
j. The Laplacian A is a bounded operator, and the spectrum of —A 
is absolutely continuous and coincides with [0,2(i]. The corresponding 
quadratic form is 

d 

(1.2) QoN = 5^5^|m(x+1,)-m(x)|2. 

Being considered on the set 5" of sequences u with finite support, the 
quadratic form Qq is non-degenerate. For n G the Hardy inequality 
is fulfilled, 

(1.3) Qo[«] > c, 5^ d>3; 

^ — ^ xr + 1 

see Section H] where we discuss this inequality and its generalizations. 
It follows from (11.31) that the completion of 5" in the metric generated 
by the quadratic form Qq is some Hilbert space of number sequences. 
We denote it by = ^{^(Z'^). The symbol V in ([H]) stands for the 
discrete potential V = V{x) > 0. We usually assume that V{x) — as 
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|x| — oo; then the operator of multiphcation by V is compact in i2{'^'^)- 
Finally, a > in fll.ll) is a large parameter (the coupling constant). 

If V —» at infinity then for any a > the essential spectrum 
aess(H.av) is the same as for a = 0, i.e., it is [0,2(i]. The negative 
spectrum consists of a finite or countable set of eigenvalues, each of a 
finite multiplicity, with the only possible accumulation point at / = 0. 
We denote by N_(liav) the number of negative eigenvalues, counted 
with their multiplicities. If V has infinite support, then A^_(Hav) oo 
as a grows. For finitely supported V one obviously has 

(1.4) iV_(H„y) < #{x G Z*^ : V(x) ^ 0}, Va > 0. 

Our main goal is to find order-sharp estimates for N_{liav)j depending 
on the properties of the potential V and on the value of a. 

A similar problem for the 'continuous' Hamiltonian is well studied, 
see e.g. the survey paper [IB] and references therein. Here we recall 
the basic results for the latter problem, since the comparison of these 
two problems will be one of our main concerns. 

The most important result for the continuous case is given by the 
Rozenblum - Lieb - Cwikel inequality. We present its most complete 
formulation, see Theorem 2.1 in [18]. This formulation goes back to 
the lectures [1], see Theorems 4.14, 4.15, and 4.17 there. For the 
Schrodinger operator on M'^ we use the same notation fll.ll) as for its 
discrete counterpart; the meaning of all terms in the formulas (II. 5p and 
(11.61) below should be clear from the context. 

Theorem 1.1. Let d > 3. Then there exists a constant C1.5 = 6*1.5(6?) 

such that for any V G (M*^), > 0, 

2 

(1.5) N_{U^v) < Ci.5«^ / vUx, 
and moreover, the Weyl asymptotic formula holds: 

(1.6) N^CHav) ~ t^dft^ / V^dx, a ^ 00. 

Conversely, suppose that d > 3, for a certain V > the operator Hav 
is well defined (via its quadratic form) and bounded below for all a > 0, 
and N^{liav) = 0{a^) as a ^ 00. Then V G Ld(]R'^), and, therefore, 
estimate (11.51) and asymptotic formula (II. 6p are fulfilled. 

We immediately conclude from Theorem 11.11 that in the continuous 
case the behavior N_(liav) = 0(0:2) is impossible, unless V = 0. The 
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growth faster than N_{liav) = 0{a2) is possible, see the results in 
[21 S] ; some of them are also presented in [IB] . 

In the discrete case an analogue of (11 .Sp remains valid, but other 
statements of Theorem II. II fail to be true. Let us formulate the result. 

Theorem 1.2. Let d > 3. Then there exists a constant Ci,7 = Ci,i{d) 
such that for any V e id{Z'^), V >0, 

(1.7) iV_(H,v')<Ci.7ai5^\/(x)i 
At the same time, 

(1.8) N_{H^v) = o{J). 

The estimate ( 11.70 is known, see e.g. [HI [IE]. The property (11.81) is 
quite elementary, but, probably, was observed for the first time in the 
survey paper [TS]. See Section [2^2] below for a more detailed discussion 
of Theorem 11.21 

Comparing these theorems, we readily see not only analogies but also 
distinctions between the continuous and the discrete cases. There are 
also others facts of a similar nature which we are going to discuss. 

Let us describe the structure of the paper. In the next Section [2] 
we present the necessary auxiliary material, and then we explain the 
proof of Theorem 11.21 In Section [3] we obtain the estimates of the 
type N_(liav) = 0{a'^) with 2q < d. They have no analogues in the 
continuous case. 

Technically, the case 2q > d is more involved than the previous one. 
The general multidimensional discrete Hardy type inequalities are an 
important ingredient of our approach. They are known in 'mathemat- 
ical folklore', but we could not find any exposition of this material in 
the literature. To fill this gap, we give such an exposition in Section 
m We describe a wide class of discrete Hardy inequalities that can be 
derived directly from their classical analogues for M'^. 

Section El is devoted to the estimates N_{liav) = 0{a'') with 2q > 
d. Based upon the Hardy inequalities established in Section HJ we 
obtain the upper estimates similar to those known for the continuous 
Schrodinger operator. We also analyze an example showing that these 
estimates are order-sharp. 

In Section [6] we analyze a special class of so-called sparse potentials. 
The estimates we derive for such potentials, depend not only on the 
values of V{x) (as, say, in (11.71) ). but also on the geometry of the 
support {x : V{x) ^ 0}. In particular, the results of this section allow 
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US to construct potentials V such that 

N^Hav) = 0{J) but N^U^v) ^ o{a^). 

By ( II. 8p . such potentials cannot lie in id. 

The results of Section[H]show also that the class of Hardy inequalities, 
described in Section HJ does not cover all the possibilities. More pre- 
cisely, we prove that any bounded sparse potential is a discrete Hardy 
weight, which may be considered as a rather unexpected result. We 
consider the material of Section [6] as the most important in the whole 
paper. 

The concluding Section [7] is devoted to a discussion of some related 
results. 

Our notation is rather standard, or it is explained in the course of 
the presentation. We denote by C,c, etc., without an index, various 
constants whose value is indifferent for us. Our notation for the more 
important constants is clear from the following example: C1.5 is the 
constant appearing for the first time in the inequality (1.5). Sometimes, 
we explicitly indicate the parameters affecting the value of C. 

The work on the paper started in April 2008, when G.R. was enjoying 
the hospitality of the Weizmann Institute of Science. G.R. acknowl- 
edges his gratitude to the Institute. 

We also express our thanks to Professor V. Maz'ya for a useful dis- 
cussion. 

2. Preliminaries. Proof of Theorem 11.21 

2.1. Birman — Schwinger operator. When studying the quantity 
N_{'Hav), it is usually convenient to move on to another, but equiva- 
lent, setting of the problem. 

Let d > 3, and let V he a bounded, real- valued function on Z'^. In 
the Hilbert space 5{^(Z'^) we consider the quadratic form 

(2.1) bvM = ^ \/(x)|m(x)P. 

If hv[u] is bounded in J{^(Z°'), it generates in this space a bounded 
linear operator, say By. We call it the Birman - Schwinger operator 
for our original problem. Since V is supposed to be real-valued, the 
operator By is self-adjoint. If > 0, it is non-negative. 

For a non- negative compact operator T, we denote by lj(T) its pos- 
itive eigenvalues, counted according to their multiplicities and num- 
bered in order of decrease. The symbol n+(s, T), s > 0, stands for the 
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distribution function of the eigenvalues: 

n+(s,T) = #{j:/,(T)>.}. 

The following statement is a particular case of the general Birman - 
Schwinger principle; see, e.g., [1] for its exposition. 

Proposition 2.1. Let V If the operator "Qy is compact, then for 
the Schrddinger operator (11.11) the number of negative eigenvalues is 
finite for any a > 0. Moreover, the equality holds: 

(2.2) A^-(H„y) = By), Va > 0. 

Conversely, if N_(liav) < oo for all a > 0, then the operator By is 
compact and, therefore, ( 12. 2p is valid. 

Recall that the modulus of an operator T is defined as 

|T| = (T*T)5. 

Evidently, |T| = T if the latter operator is non-negative. The compact 
operators, such that for some g > 

n^{s,\T\) = 0{s~^), 

or, equivalently, /j(|T|) = 0(j~^), form a linear space (an ideal in the 
algebra of all bounded operators). This space is usually denoted by 
&q,oo, or by Eg. The latter notation is used in the book [2], and we 
keep it here. The spaces are the 'weak analogues' of the classical 
Schatten ideals &q. The functional 

(2.3) ||T||s, = sup sn+is, |T|)^ = sup j'/,(|T|) 

s>0 j 

defines a quasi- norm on S^. If g > 1, and only in this case, an equivalent 
norm exists on Eg. However, the quasi- norm (12. 3p is more convenient 
for estimates. 

The spaces Eg are complete and non-separable. The condition 
n+(s,|T|) = o(s-'?), 8^0, 

or, equivalently, /j(|T|) = o(j~^), singles out a closed separable sub- 
space of Eg, which we denote by E°. See Section 11.6 in the book [2] 
for more detail on the classes Eg and E°. 

Along with the quasi- norm (12.31) . let us consider the functionals 

Ag(T) = limsup(s%+(s, T)), 5g(T) = liminf (s%+(s, T)). 

They are well-defined for any T G Eg. The following elementary fact, 
see Theorem 11.6.7 in [2], will be useful later. 
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Proposition 2.2. The functionals Ag, 5q are continuous in the topology 
of the space S^. 

Proposition 12 . 1 1 allows one to replace the study of the function N_{\iay) 
with the investigation of the compact operator By. This is convenient, 
since it makes it possible to use the powerful machinery of the theory 
of compact operators. 

For instance, for any g > the powerlike estimate 

iV-(H«y) < Ca\ Va > 

- 

is equivalent to By G Sg, with ||By||2 <Ci. In a similar way. 



9' 

iV_(H,V') =oK) ^ ByGS°. 



In particular, the estimate fll.7l) . being reformulated in terms of the 
operator By, reads in either of two equivalent forms: 

(2.4) ||By||s, < CI,\\V\\,^- n+(s,By) < C^.-js-i V V{x)i. 

T. T. 

The spaces Eg and S° have their counterparts Iq^^ and £° in the 
theory of function spaces on countable sets. Let a function V on I/- 
be such that V{x) as \x\ — > oo. By re- arranging the numbers 
X G Z'^, in the non- increasing order, we obtain a sequence of 
non- negative numbers (notation V* , j G N). Then V^* — as j — > oo. 

Besides, for any s > we denote 

(2.5) E{s, V) = {xeZ'^: \V{x)\ > s}; u{s, V) = i^E{s, V). 
By definition, 

{veiq,^} ^ {u{s,v) = ois"')} ^ {v; = o{r')}, 

and 

{veeu ^{u{s,v) = o{s-^)} ^ {v; = oirl)}. 

The functional 

ll^ll^g^ = sup su{s,V)^ = sup j^V* 

s>0 j 

defines the standard quasi- norm in ig^^. 

The next statement elucidates the role of classes ig^^ and £° in the 
problems we are studying. For its proof, see Section 14.21 below. 

Proposition 2.3. Let d > 3 and V G ^|,^(Z'^), 1^ > 0. Then the 
operator By is bounded. IfVEl\ (Z"^); this operator is compact. 
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Note that the conditions given by this proposition are only sufficient, 
but not necessary. Say, they are violated for the sparse potentials 
considered in Section [61 

2.2. On Theorem 11.21 As it was already mentioned in Section [H 
Theorem 11.21 is basically known. Still, in order to make our exposition 
self-contained, we outline its proof. 

Proof of Theorem It is well known that the discrete Laplacian 
generates a positivity preserving semigroup in the space If 
(i > 3, the following lower estimate, with some c > 0, is satisfied for 
the quadratic form Qq in (11.21) : 



This is a particular case of Theorem IV. 5. 2 in [T^. Thus, the assump- 
tions of Theorem 1.2 in [TT] are fulffiled, whence the estimate (11. 7p . 

Another way to obtain the latter estimate is to derive it from the 
results of the paper [17] ; see also Section 9 in [18] . It is also possible to 
deduce (I1.7P directly from (11.51) . This will be explained in Section [5TT1 

For justifying (11.81) . let us consider the operator By. The sequences 

V with finite support form a dense subset of £<i(Z°'). By (II. 4p . for such 

V the non-zero spectrum of By is finite and therefore, ( 11. 8p is satisfied, 
so that Ad (By) = 0. Using the inequality (12.41) and Proposition [221 we 

conclude that (II. 8p extends to all non-negative potentials V G idiJL^). 
The proof of Theorem 11.21 is complete. 

Note that the estimate (II. 8p does not replace the, formally weaker, 
estimate (11.70 . An important property of (11.70 is the explicit depen- 
dence on the function V . It is this dependence that enabled us, by 
approximating V by potentials with finite support, to pass to the limit 
and obtain (II. 8p . 



As it is clear from the comparison of Theorems 11.11 and 11.21 the esti- 
mates obtained in this section have no analogues for the operators on 
M"^. We are going to show that, unlike in the continuous case, the func- 
tion A^„(HQ,y) can grow as 0{ofl^ with arbitrarily small g > 0, which 
corresponds to the arbitrarily fast powerlike decay of the eigenvalues 
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Theorem 3.1. Let d > 3, q < f , and V G ig,yj{Z'^), V > 0. Then 
By G Eg, and 



(3-1) l|Bv||s, < C3.1I1 V 

with a constant depending only on d and q. Equivalently, 

N_{Yi^y)<Cl,anV\\l^. 

For = 1 the statement is no more true: the inclusion V E id , 
guarantees only the boundedness of By. 

Proof of Theorem \3.1[ For a fixed s, let us split V into two terms, 

V = Vi + V2, 

where Vi(x) = V{x) at the vertices x G E{s, V), cf. (12.51) . and Vi(x) = 
otherwise. By the Weyl inequality (see, e.g.. Theorem 9.2.9 in |2]), 

(3.2) n+(s. By) < n+{s/2, ByJ + n+{s/2, ByJ. 
The first term in (13. 2p is estimated by (11.41) . 

(3.3) n+(s/2, ByJ < z/(s/2, V) < {2P/s)\ P = 



For the second term we apply the estimate (12.41) . with V replaced by 
V2. This gives 

(3.4) n+(s/2,ByJ<2fCi.7s"f Yl ^(^)^- 

V{x)<s 

Now we transform the latter sum: 

J2 r(x)i = - rrldu{T,V) = ^ r T-^'^iy{T,V)dT 

V(x)<s ^ 

- 2 Jo - 2g 

Together with (D, (O, and ([32D, this gives ri(s. By) < CP^s"" 
which is equivalent to (13.11) . □ 

It is possible also to obtain the lower estimate of r;,+ (s. By) in terms 
of the function z/(-, V). This estimate does not require any preliminary 
assumptions about the non-negative function V . 

Theorem 3.2. For any V >Q the inequality holds: 

(3.5) n+(s. By) > 2-'^u{2sd, V). 
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Proof. For a fixed r, we consider the set E{t.,V). Let the sublattice 
{2'LY consist of the points in with all components divisible by 2. 
For any k G {0, l}*^ consider the set (2Z)°' + k. These sets are mutually 
disjoint. Therefore, the sets f2k(V^) = E{t, V) fl ((2Z)'^ + k) are disjoint 
as well, and at least for one value ofkG{0,l}'^, we have 

i^n^{y)>2-''v{T,v). 

Now, consider the subspace L C !K^(Z"') formed by the functions 

u{x)= ^ Cy6{x-y). 

The functions 5{x — y), y E Q\t_{V), are mutually orthogonal both in 
the metric of and with respect to the quadratic form hy in 

fl2.1l) . So, for any u G XL we have Qo[u] = '2dJ2yen^{v) l^yl"^^ while 
byM = Eyefik(v) IcyP^d/)- Since V{y) > r for y G ^k{V), we have 
constructed a subspace of dimension greater than 2~'^i/(r, V) such that 
by['u] > {2d)~^TQo[u]. This immediately implies (13. 5p by the varia- 
tional principle. □ 

Corollary 3.3. Let d > 3, < q < ^, and V > 0. Then By G if 
and only ifVE Iq^^ilJ^), and, moreover, there are constants Co,Ci > 
such that 

(^oW\\tq,Ai<^) < llBylls, < Ci||V^||£,,„(zd), 

and also 

(3.6) Co limsup s^i/(s, V^) < Aq(Bv) < ci limsup s''z/(s, V^). 

In particular. 

Here only the relation (13. 6p needs a justification. To this end, let us 
denote ^ 

R{y) = limsup sz/(s, V)'^ , 

then R{V) < \\V\\i^^. For any e > 0, one can change the values of 
V{x) at a finite number of points a; G Z"' in such a way that for the 
new potential, say V^, we have 

mie,^^, < RiV) + e. 

Then AqiByJ = Ag(By) and, by dSH), 

A,{By)<\\ByJ\l^<Cl,{R{V)+er. 

Since e is arbitrary, we come to the right inequality in (13. 6p . The left 
inequality is a direct consequence of (13. 5p . 
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4. Discrete Hardy inequalities 

Here we collect the material on the discrete multidimensional Hardy 
type inequalities. We need them for studying the estimates of the op- 
erator By in the classes with q > ^. As it was mentioned in the 
Introduction, this material should be considered as known on the 'folk- 
lore level'. However, we could not find its exposition in the literature, 
and decided to present it here. 

In this section we describe a rather simple class of discrete Hardy 
inequalities. Namely, with any function W{x) > on Z"^ we associate 
a function W(^) on M"^, which assumes the value W{x) on the unit 
cell determined by the vertex x; see Subsection 14.11 for details. Then 
we show that if the "continuous" Hardy inequality is satisfied with 
the weight W(^), then the discrete Hardy inequality with the weight 
W{x) holds on Z*^. Naturally, this class of discrete Hardy inequalities 
is rather restricted. Later, in Section 16.31 we will show that there exist 
discrete Hardy inequalities of a different origin. 

4.1. Operator of poly-linear interpolation. Let 'K^{W^),d > 3, 
stand for the homogeneous Sobolev space. The metric in !K^(M°') is 
defined by the standard Dirichlet integral, 

D[u] = I iv^(orrfe. 

We introduce an operator J : J{^(Z'^) 'K^{W^) in the following 
way. In every elementary cubic cell in R'' we interpolate 2'^ values of u 
at the vertices of the cell by a poly-linear function. For instance, the 
function u{x)^ x G {0, 1}'', such that m(0, . . . , 0) = 1 and u{x) = at 
the remaining vertices, interpolates as 

f/(0 = = - 0), i = (ei, [0, i]'^. 

i 

Such interpolation defines a mapping : u ^ U of the space "K^ilJ^) 
to a certain subspace of piecewise poly-linear, continuous functions on 
W^. 

Lemma 4.1. The image of the mapping J is the space 'Ky^{W^) of all 
piecewise poly-linear, continuous functions with the finite Dirichlet in- 
tegral. The quadratic forms Qo[u], see { \1.2\i . and DlJu] = f^^ \V{'Ju)\'^dx 
are equivalent: there are constants c, c' > such that 

(4. 1) cQo [u] < D pu] < cQo [u] , Vm G (Z'^) . 
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Proof. Consider the space Cj{C) of poly-linear functions on the unit cell 
C = [0, 1]"'. Clearly, dim £(6) = 2^. On £(C) we consider the quadratic 
forms 

Q[U; e]= Yl - uiy)\'; d[U; e] = I |Vf/(Oprfe 

These two quadratic forms vanish on the same subspace in 'C(C), con- 
sisting of constant functions. Therefore, they are equivalent, i.e., with 
some c, c' > we have 

cQ[U; e] < D[U;e] < c'Q[U;e]. 

By adding up similar inequalities for all the cells G + x, x E Z*^, we 
arrive at (14.11) . 

It remains to check that Ju G J{^(M"') for any u G J{^(Z°'). It is 
sufficient to show this for the dense in "K^ [7/) subset of functions with 
finite support. But for any such u the function has compact support 
and, therefore, it can be approximated in "K^ by functions in C^, in 
the metric of the Dirichlet integral. Hence, it lies in J{^(M'^). □ 



Now, with any non-negative function W{x), x G Z"', we associate a 
function W = on M*^, setting W(,^) = W{x) for ^ G Q + x where C 
is the same as above. 

Lemma 4.2. For any weight W on and the corresponding weight 
W = 3W on R'^, one has 

(4.2) 5^H^(x)Kx)p<C4.2 / W{0\U{0\'dt U = 3u, 

for any function u onU^ . The constant in (14.21) depends only on d. 

Proof. By linearity, it is sufficient to prove (14. 2 p for W{x) having sup- 
port at one point, say, at G Z'^. Then the quadratic forms in (14.21) 
define Hilbert seminorms on a 2'^-dimensional space, with the null sub- 
space for the quadratic form on the right-hand side contained in the 
null subspace for the one on the left-hand side. This gives the required 
inequahty. □ 

Note that the estimate inverse to (14. 2 p is impossible. Indeed, the 
null subspace for the quadratic form on the left-hand side is strictly 
wider than for the other one. 
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4.2. Hardy type inequalities. We will call a non-negative function 
W{x), X G Z'^ a Hardy weight on Z'^, if for some constant H{W) the 
inequality is satisfied: 

(4.3) ^ iy(x)|u(a;)|=^ < if(iy)goM, Vu G ?{^(Z'^). 

The best possible constant H{W) in (14.31) will be called the Hardy 
constant for W. We will say that a Hardy weight W is normalized, if 
if(W^) = 1. This definition carries over to the lattice case the classical 
definition of Hardy weights in M*^: 

[ W{^)\U{0\^d^ <H(W) [ \VU\^d^, WUe^\R'^). 

In the continuous case the complete description of Hardy weights was 
found by Maz'ya, see [12], Ch. 8. There, the necessary and sufficient 
condition on a function W(,^) to be a Hardy weight in M'^ is given in the 
terms of the capacity. By means of the well-known relation between 
capacity and measure, a sufficient condition for W to be a Hardy weight 
can be expressed in more elementary terms, see, e.g.. Proposition 5.1 
in [4]. 

Proposition 4.3. Any function W > in the weak class Ld ^(M*^), 
d > 3, is a Hardy weight on with Hardy constant if(W) satisfying 

HiW)<C{d)\\'^\\L, (M^)- 

The discrete analogue of Proposition 14.31 also holds. 

Proposition 4.4. Let W > be a function defined on Z'^, d > 3. 
Suppose that W = 3W is a Hardy weight on M'^. Then W is a Hardy 
weight on Z^, and H{W) < CH(W). In particular, any W G id ^(Z'^) 

is a Hardy weight on with H{W) < C{d)\\W\\i (id). 

Proof. The first statement follows immediately from the inequalities 
(14.11) and (14. 2p . As for the second statement, note that the functions 
W on Z'^ and W = on are equimeasurable, therefore W G 
Ld ^(M*^), with the same quasi-norm, so we can use Proposition 14. 3[ □ 

The simplest example of a function in id ^{Z'^) isW{x) = {\x\'^ + l)~^. 

By Proposition 14.41 it is a Hardy weight on Z'^. This justifies the 
inequality (11.31) . 

Now we are in a position to prove Proposition 12. 3[ Indeed, the 
boundedness of the operator By is just a re-formulation of the property 
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of \^ to be a Hardy weight. Evidently, By is compact for any V with 
finite support. This property extends to the whole of by continuity. 

2,W 



5. Estimates in S^, g > f 

Having the Hardy type inequalities at our disposal, we now move 
on to the estimation of the operator By in the classes Eg with q > ^■ 
Results in this section are the direct analogues (actually, immediate 
consequences) of the corresponding results for the operators on R"^, 
obtained in [3lll]. In order to distinguish between the discrete and the 
continuous Laplacians, we (in this section only) will denote them by 
Arf and by Ac respectively. 



5.1. Eigenvalue estimates. The material of Section H] allows one to 
prove that any statement on the eigenvalue behavior of the operator 
By, expressed in terms of the standard, or weak Lg-classes, automati- 
cally implies its discrete counterpart. In particular, this shows that the 
estimate (11.71) follows directly from (11.51) . thus giving one more proof 
of Theorem II. 2[ Below we demonstrate, how this idea implements for 
obtaining the discrete analogues of the results in [3l 0] . 

Let a function W be a normalized Hardy weight on M°', d > 3. The 
following is the formulation of Theorem 4.1 in [1]. For a function F > 
on we say that F e Lg,^(W2) if 



|Ff , =sup 

Lg.n,{^^) t>0 



( \ 

\ F(C)>tW{0 / 



< oo. 



Proposition 5.1. Let d > 3, and let W > be a normalized Hardy 
weight on W^. Suppose | < g < oo and V > 0. Then 

V/WgL,,^(W'^/2) ^ iV_(-A,-aV)<C5.i«^||V/W||^^_^(^,/,^ 

where the constant depends on d and q. 



If, besides, Jim | t« / W^/'^d^ ) = 0, t/ien iV_(-Ac-aV) 

V(0>tW{0 / 



We now can prove the discrete version of this Proposition. Let W > Q 
be a Hardy weight on Z*^, d>3. For a function F > on TJ^, we say 
that F e ^q,^{W'^/^) if 
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\ F(x)>tW(x) 

Theorem 5.2. Let d> 3, and letW > be a function on , such that 
W = is a normalized Hardy weight on W^. Suppose | < g < oo 
and > 0. Then 
(5.1) 

M/ c P (U/'i/'^\ ^ AT f—A . ^ <r r. .o,9|lT//U/l, 



Ij, besides, lim U« ^ ^^(x)^/^ | _ g, t/ien A^_(-Ad-a1/) 

t^0,oo \ y(^)>jH/{a;) 



Proof. Consider the operator Bv in the space IK^(]R'^), defined by the 
quadratic form 

bv[f/] = j v(oit/(OI'rfe, v = av^. 

By the Birman-Schwinger principle, see Proposition 12.11 (or, more ex- 
actly, its continuous version), we have 

- aV) = n+{a-\ Bv). 

The function V is equimeasurable with V, W is equimeasurable with 
W, and V/W is equimeasurable with V/W. Therefore the condition 
V/W e ig,^iW^/^) implies V/W e L,,^(W'^/2), so we can apply The- 
orem 15.11 and obtain the estimate 

(5.2) n4a-\ Bv) < C,.,a^V/Wh^^^^^.f.y 

Now we consider the operator Bv defined by the same quadratic form 
bv[?7], but restricted to the space 'K7^(M.'^) of piecewise poly-linear 
functions in J{^(M'^). Since we narrowed the domain of the quadratic 
form, the eigenvalues of the operator cannot grow, so fl5.2p leads to 

(5.3) t2+(«-i,B0.) < Cs.iall Wllv^CH--/^)- 

By (14.1 p and (14.20 . the quadratic form hy[u] in (12. ip is estimated 
from above by bvpw], and the quadratic form Qo[u] in (11.20 is esti- 
mated from below by J \'V{du)\'^dC,. Therefore the eigenvalues of the 
operator By are majorized by the eigenvalues of the operator Bv, and 
the estimate (15. ip follows from (15. 3p . The second part of Theorem 15.21 
is proved by a similar reasoning. □ 
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Theorem 15.21 gives much freedom in choosing the Hardy weight W. 
The most standard choice is W{x) = c(|xp + l)^""^ where c > is the 
normahzing constant. For this W, the corresponding particular case of 
Theorem 15.21 was presented in [18] as Theorem 9.2. Below we repeat 
its formulation, in order to correct a misprint in [TSj . 

Corollary 5.3. Let d > 3 and 2q > d, and letV Then 

(5.4) Ar_(H„v) < C(d,g)a«sup I ^ (Ixp + 

\ (\x\-^+l)V(x)>t 

5.2. An example. Here we analyze an example which shows that the 
estimate (15. 4p is order-sharp for certain potentials V . Note that in the 
setting of Section [3] such examples are unnecessary, since its results give 
two-sided estimates of the quasi-norm ||By||sg for 2q < d. 

What we give below, is an analogue of Example 4.1 in pj. However, 
in the discrete case we were not able to establish the asymptotic be- 
havior of the eigenvalues lj(Bv), and we only give for them some lower 
estimate. 

In our example the potential, for |x| > 1, has the form 

V{x) = \x\-^{\og{\x\))~^, 2q > d. 

By Corollary 15. 3[ for the corresponding Birman-Schwinger operator 
By, the upper eigenvalue estimate holds, 

n+(s,By) < 

We will show that this estimate cannot be improved. More specifi- 
cally, we will prove that By ^ S°. To this end, we will construct an 
orthonormal in sequence of functions (pn such that 

(5.5) by[(y9„] > era" 9, c > 0. 

Then the desired lower estimate will follow from the well known re- 
sults on compact operators in a Hilbert space. Namely, it follows from 
Lemma II. 4.1 in |8], or Theorem 11.5.7 in [2], that 

m m 

$^/„(By) > ^by[(^„], VmeN. 

n=l n=l 

Since ^ < | < 1, this inequality and (15.51) imply 

m 

^/n(By) > Cm^"9, 
n=l 

which is inconsistent with By being in the class S°. 
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So, it remains to construct the functions <y9„. We obtain them by nor- 
mahzing a certain sequence {un}- For constructing the latter, we use 
the Fourier representation of functions on the lattice. With a function 
u G "K^ we associate its Fourier transform, 

u{z) = (27r)-'^/2 J2 M(x)e-^^■^ z G 
where T*^ stands for the (i-dimensional torus. Conversely, 



uix) = (27r)-^/2 / ^z)e''''dz. 



Then 



~ d 

(5.6) Qq[u]= \u{z)\'^uj{z)dz, uj{z) = 4:S2sm^{zj/2) 

.,=1 



We take = 4 " and set 

.„W = s..W = x-(g-3)nx(^). 

where x is the characteristic function of the interval (—1, 1). So, the 
function Vn{z) has support in the square with the side length 2/i„ and 
with the center moved by 3/i„ in the direction of zi axis. Obviously, 
the supports are disjoint, so the functions f„ are orthogonal in L2,, 
The latter notation stands for the L2-space with the weight uj. Thus 
the functions m„ are mutually orthogonal in "K^. 

The functions m„ can be calculated explicitly: 



5 



(5.7) Un{x) = (27r)-'^/2g3^h„xi -Q 



■ 1 '^i 



We have ||ttn||j<;i = C*|bn||i2„- Since uj{z) x h'^ on the support of m„, 
this gives ||w„||i,^ X hlWvJl^ X hl+'^. 
Now, from (15. 7p we derive that 



We need to estimate the last expression from below. To this end, we 
restrict summation in the last sum to the region f2„ C Z'^ where all 
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\hnXj\ < 7t/2 and |x| > 4. In fi„, we have sin^(/i„Xj) > (^^hnXj^- 
Therefore, 

hv[Un] > Chl' ^(^) = c/^n' E l^r'^^S I^D^ 

A lower bound is given by the integral 

rch-^ 1 1 



^4 

Therefore, 

by[Mn]/||Mn||^i >c|log/i„|~9 >cn"9. 
So, by normalizing the functions Un we obtain the sequence {v^n} sat- 
isfying (15. Sp . 

It is interesting to notice that the test functions guaranteeing the 
lower estimate for 2q < d were constructed as having disjoint supports, 
while for 2q > d such functions have disjoint supports of their Fourier 
transform. 

6. Sparse potentials 

The results of Sections [3l [5] allow one to construct, for any prescribed 
value of g 7^ |, the potentials V such that By G Sg, but By ^ E°. For 
the borderline value q = ^, Theorem 11.21 leaves open the question of 
existence of such potentials. 

In this section we consider a special class of potentials for which the 
theory can be advanced much further. In particular, we answer the 
above question by showing that for any number sequence pj \ such 
that pj+i/pj ^ 1, a potential V does exist, such that the sequence of 
eigenvalues lj(Bv) asymptotically behaves as {pj}- 

6.1. Green function of the discrete Laplacian. The operator (—A)^ 
acts as a discrete convolution, its kernel can be represented by the ex- 
plicit formula: 

hy{x) = ho{x - y), \/y E Z'^, 

where 



(6.1) /.o(x) = (2vr)-M , . dz. 

Jt'' 4 2]j=ism {zj/2) 

Note that here the denominator coincides with the weight function u 
in (ra . 

We will call ho the Green function. One should be careful when 
studying its properties, since the point I = lies in the spectrum of 
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—A. However, all the difficulties can be easily overcome by systematic 
use of the representation fl6.ip . See, in particular, [5j, where the case 
c? = 3 is analyzed. For any d > 3 the reasoning is similar. 

The function hQ{x) lies in [K^(Z'^), is harmonic outside the point 
a; = 0, and its value at this point is 

M0)=/.^ = (2vr)- / , ■ 

Then also hy{y) = /io(0) = /i^ for any y G Z"^. Besides, 

{Ahy)iy) = -1, yyeZ". 

Let M be a function with finite support. Then summation by parts 
leads to the equality 

(6.2) (W, hy) = -Y. n{^){^hy)i^) = U{y). 

Starting from (16. 2p . the scalar products, and also the norms, are taken 
in "K^ . Since finitely supported functions are dense in CK^, the equality 
( 16.21) extends by continuity to all m G In particular, 

(6.3) {hy,hy^) = hy{yi) = ho{y ~yi). 
Taking in (16. 3p yi = y, we find that 

\\hy\\ = fi, G Z'^. 

We have < hQ{x) < for all x G Z'^. It also follows from the 
representation (16.11) that 

(6.4) ho{x) < CeA\x\-^^-^\ x ^ 0, 

with a constant depending only on d. The relations (16. 3p and (16.41) 
show that for the points y, yi G Z'^ lying far enough from each other, the 
functions hy, hy^ are 'almost orthogonal'. It is convenient to normalize 
these functions, so that further on we work with 

6.2. Sparse subsets in Z"'. Let y be a subset in Z*^, d > 3, and 
let Oiy stand for the subspace in 'K^{l/') spanned by the functions 
hy, y eY. We say that the set Y is sparse (or, in more detail, strongly 
sparse), if in 'K^ilJ^) there exists a compact operator T, such that the 
operator I — T has bounded inverse and the functions 

(6.5) ey = {l-T)-%, yeY, 
form an orthonormal system in IK^. 
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We say that Y is weakly sparse, if in the above definition we replace 
the requirement of T being compact by its boundedness. 

We are going to describe (in Lemma ESI below) a rather general way 
to construct sparse subsets. To this end, we need some preliminary 
material. 

Suppose that a sequence {'i/'j}, j G N, of elements of a separable 
Hilbert space can be represented as 

where {^j} is an orthonormal basis in and T is a linear operator of 
the class ©g with some g, < g < cxd. Then we say that {ipj} is a 
g-basis. We recall that ©oo stands for the space of all compact linear 
operators in S), and &q, < g < oo, stands for the Schatten ideal; see, 

e.g., mn- 

The following result is due to Prigorskii, see Theorem 5 in [16]. 

Proposition 6.1. Let {ipj}-, j & N, be a complete and uj-linearly in- 
dependent sequence in S). It is a q-basis if and only if the matrix 

belongs to the class &q. 

Here the property of cu-linear independence of the system {ipj} means 
that the assumptions 

'n = {vj}^i2, ^Vj'^j = o 

j 

yield rj = 0, see definition in [8], Section VI. 2. 4. 

Note that for g = 2 the statement of Proposition 16.11 turns into that 
of Theorem VI.3.3 in [8]. The proof for the general case follows the 
same scheme as in [8]. 

Now we introduce the quantities that appear in the formulation of 
Lemma [6.21 below. Let y be a subset of Z*^. Given a. y eY, we denote 

Ty = dist{y,Y\{y}) 

and 

[y] = #{x G F : |x| < \y\}. 

Lemma 6.2. Suppose d > 3, and let Y d be a set, such that 
[?/]rj7^'^~^'' — > as ?/ G y, \y\ ^ oo. Suppose also that 

J2 ^ ^'"'^ < ^' sup {[y]r-^'-'^) <A, Vy G Y, 
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with a constant A satisfying 
(6.6) 2AC6A < /ul 

Then the set Y is strongly sparse. 
Proof. Consider the Gram matrix 

G = ((h^Ji 



x,y£Y 



x,yeY; \x\,\y\<N 



Its diagonal elements are equal to one, and by (16.31) . for every y E Y 
we have 

\ihx,hy)\< C6.4^-M E + E ) |a;-y|-('^-2) 

x&,xjty \^(^y\{y},\x\<\y\ xeY,\x\>\y\J 

(6.7) < Ce.,f^-'([y]ry^'-'^+ E r-^''"^ 

y x£Y,\x\>\y\ 

Since the matrix G is Hermitian, the same inequality holds if the 
roles of X and y are interchanged. By the Shur test (see, e.g.. Theorem 
2.5.6 in [2]), this implies that the matrix G — I defines a bounded 
operator in l2iX) and, moreover, (16.61) yields \\G — I\\ < 1. 

In order to prove that G — I E ©oo, we show that the 'truncated' 
matrices 

Gn = {ih.j., hy 

converge to G in the operator norm. For estimating HG — Gjvll we 
again use the Shur test. For the rows with \y\> N the estimate (16.71) 
survives, and under the assumptions of Lemma the sums on the right 
are small if is large enough. If \y\ < N, then this sum should be 

replaced by E ^'^ '^^ which is also small. So, G — I E &oo- 

x&Y,\x\>N 

Next, we check that the system {hy]y^Y is ci;-linearly independent. 
Indeed, by taking the scalar products with /i^, x G F, we come to 
the infinite system of linear homogeneous equations, Gr] = 0. Since 
||G — J|| < 1, the operator G is invertible in £2, so that the only £2- 
solution of this system is trivial: 77 = 0, and we are done. 

Now Proposition 16.11 applies, with q = 00. So, we get an operator 
To G ©00 in the space "Ky, such that the system {(I — To)/iy}j,ey 
is orthonormal. Extending Tq by zero to the orthogonal complement 
of "Ky, we obtain the operator T that satisfies all the properties we 
need. □ 

We do not have a method that would allow us to construct subsets 
that are weakly sparse but not sparse. Still, we consider the notion 
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of weak sparseness useful, since it gives us a way to see the difference 
between implications of both types of sparseness. 

6.3. Sparse potentials. We say that ^ > is a sparse (a weakly 
sparse) potential on Z"', if its support 

Yv := {xeZ'^: V{x) > 0} 

is a sparse (respectively, weakly sparse) subset. 

Let be a weakly sparse potential. Then, using the equalities 06.21) 
and (16.51) . we see that the quadratic form hv[u] can be written as 

hy[u] = J2 hy)\' = f.'Yl ^(^)i(^' (I - T)e.)r- 

y^Yv y&v 

Along with the operator By generated by this quadratic form, consider 
also the self-adjoint in "K^iT/) operator 

(6.8) ^y = ^Y, Vny){-.ey)ey. 

y&Yv 

We have 

Ny(I-T*)=/i5^ V^(-,Me„ 

whence 

||Ny(I - T*)wf = /i' 5^ V{y)\{urhy)? = hy[u]. 

y&v 

This means that 

(6.9) By = (I-T)N^(I-T*). 

This representation allows one to obtain the following simple result 
that is not covered by Proposition 14.41 

Theorem 6.3. Let V > Q he a function on Z'^, d > such that its 
support Yy is weakly sparse. Then V is a discrete Hardy weight if and 
only if V is hounded. Moreover, the following two-sided inequality is 
satisfied for the Hardy constant H{V): 

(6.10) C\\Vl^<HiV)<C'\\Vh^. 

where C = ||(I - T)-i||-V and C = ||I - T|| V- 

Proof. In the equivalent terms, we have to find the boundedness con- 
ditions of the operator By. They immediately follow from the repre- 
sentation (16. 9p and the equality ||Ny|| = . □ 
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The next result is also a consequence of (16.91) . Recall that for a 
function V, such that V{x) as |x| oo, we write V*,j G N, for 
the numbers |V^(a;)| rearranged in the non-increasing order, see Section 

Theorem 6.4. Let V > be a function on Z*^, d > 3, such that 
its support Yy is weakly sparse. Then the operator By is compact if 
and only ifV{x) ^ as \x\ oo. Moreover, the following two-sided 
inequality is satisfied for the eigenvalues lj(By): 

CV* < I J (By) < C'V*, Vj G N, 

with the same constants as in (16.101) . 

The proof is the same as for Theorem 16. 3[ 

The next result is more advanced, and it requires the potential to 
be sparse but not weakly sparse. We derive this result for potentials V 
which meet an additional condition: we assume that the correspond- 
ing sequence {V^*} is moderately varying. We say that an (infinite) 
sequence {pj} of positive numbers is moderately varying, if pj \ and 
Pj-^i/pj 1. We shall use a result of M.G. Krein, see Theorem 5.11.3 
in [8]. Below we reproduce its formulation, restricting ourselves to the 
situation we need. 

Proposition 6.5. Let H > and S be self-adjoint, compact operators. 

Suppose rankH = oo and the sequence of non-zero eigenvalues 

is moderately varying. Then for the operator M = H(I + S)H one has 



lim -4r, — r = 1- 
/|(H) 

The next theorem is the main result of this section. 



Theorem 6.6. Let V > be a sparse potential, such that the numbers 
V* form a moderately varying sequence. Then 

lim MM) = 1. 

Proof. The non-zero spectrum of the operator (16. 9p coincides with that 
of the operator 

(6.11) My := Ny(I - T*)(I - T)Ny = Ny(I + S)Ny, 

where 
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Now we apply Proposition 16.51 to the operators H = and My, 
given by fl6.8l) and fl6.1ip respectively. All the assumptions of proposi- 
tion are evidently satisfied, and we get the desired result. □ 

Corollary 6.7. For any moderately varying sequence {pj} there exists 
a sparse potential V > 0, such that lj(Bv) ~ Pj- 

Taking here pj = with an arbitrary g > 0, we obtain a potential 
V such that the eigenvalues /j(By) behave as pj and, therefore. 

For g = |, this solves the problem stated in the end of Introduction. 
This corollary allows one also to construct potentials V , such that the 
eigenvalues /j(Bv') decay arbitrarily slowly, say, /j(By) ~ (logj)""^. 
This corresponds to the exponential growth of the function N_{H.av)- 
Applying Theorem 16. 4^ we easily obtain also potentials with the loga- 
rithmic growth of this function. 

7. Discussion and concluding remarks 

7.1. Let us briefiy summarize our main results. 

1. Each estimate of the type 

(7.1) Ar_(H,y) = 0(a«) 

with 5' > |, known for the continuous Hamiltonian (11. ip . has an ana- 
logue for its discrete counterpart. What is more, such analogue can be 
derived directly from its prototype. 

2. For (? > ^ there are also estimates (17. ip of a different origin, that 
hold for the discrete Hamiltonian (II. ip but do not have continuous 
prototypes. These are estimates for the sparse potentials. Unlike the 
results mentioned above, these ones give also lower estimates of the 
same order. 

3. Estimates of the type (17. ip with g < |, that hold for the discrete 
Hamiltonian, have no continuous analogues either. They give necessary 
and sufficient condition on the potential, in order that By € Sg with 
a given q to be satisfied. 

7.2. Here we discuss the case of non-sign-definite potentials V ^ (). 
As usual, we denote 

V^ = \{\V\±V). 

For the operator 'Ray with such potential, the eigenvalues may appear 
also to the right of the point 2d = max(Tess(HQ,v), and along with 
N-iiiav) one should consider the function N+{2d; Hav) defined as the 
number of such eigenvalues, counted with multiplicities. 
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Consider the mapping 

r : uix) ^ (-l)"^+-+"''u(a;). 
It defines a unitary operator in £2(^'^), and it is easy to see that 
r(-A - aV) = (A + 2dl - aV)T, 

therefore, 

(7.2) Ar+(2(i;H„v) = iV-(H_,v.). 
The inequahties 

(7.3) iV_(H,y) < iV_(H,^J, iV+(2(i;H„y) < iV.(H,yJ 

reduce the problem of estimating the functions on the left to the similar 
problem for the non-negative V. The first inequality in (17. 3p follows 
directly from the variational principle. To derive the second, one uses 
also the equality (17.21) . 

Sometimes, lower estimates of the same type are also possible. In 
particular, this is the case for the results in Section [31 

7.3. Up to now, the problem of counting the number of negative eigen- 
values of the multidimensional discrete Schrodinger operator did not 
attract much attention of the specialists. The only paper we are aware 
of, is [ID], where the operator (1.1) with a = 1 was considered, under 
the assumptions that can be formulated as follows: 

V = V,- V2, V^, V2 > 0, 

and, moreover, 

V2>Uo> 0; Vi eir, r> 1. 
The main result consists in the estimate 

(7.4) iV_(Hv)<(e/z/o)^5^V^i(a;)^ 

that holds for all c? > 1, rather than only for d > 3. In the latter case, 
the estimate fl7.4p is much weaker than any of our estimates. 

Below we analyze the estimate (17.41) . The conditions on V imply, 
for each n G Z'^, that V{n) < Vi{n) — uq. Hence, V^{n) > yields 
Vi{n) > vq. The number !N of such indexes n is always finite, and 
moreover. 

Using now (17.31) and the trivial inequality (11.41) , we come to an estimate 
similar to (17.41) . but without the redundant factor e. 



SCHRODINGER OPERATOR ON 



25 



In other words, the estimate (17.41] needs nothing for its vahdity, ex- 
cept for the non-negativity of the operator —A. In particular, this 
explains why it holds in any dimension. 

7.4. A problem, closely related to the one of counting eigenvalues, 
concerns the Lieb - Thirring type inequalities. This problem was con- 
sidered in [9]. Let E~{iiv), Ej~{liv) stand for the eigenvalues of the 
operator Hy, lying on (— oo, 0) and on {2d, oo) respectively. The fol- 
lowing estimates were established in [9J for any d > 1, along with other 
interesting results: if G dq^'L'^) for some (? > 1, then 

\Ej{^v)\'' + Y.^Ej{Yiy) - 2dy < J2 

j j xGl,'' 

if \^ G i^^d{Z'^), again for some g > 1, then 

J2 lEiiHv)]'^ + J2(^l{liv) - 2dr < C{d, q) J2 \v{x)rK 

with some, explicitly given constant factor. 

In this respect we note that from our results it easily follows that in 
dimensions > 3 a similar inequality holds for any g > 0, but with an 
additional and unspecified factor C'{d,q) on the right. Moreover, in 
the second estimate the assumption V G ig^d{Z^) can be replaced by 

V G ig+r{'^'^) with an arbitrary r G (0, |], and with the corresponding 
change in the right-hand side. 

In this paper we do not look for the optimal constants in such es- 
timates. Note that traditionally just seeking the optimal constants is 
the main issue when dealing with the Lieb - Thirring inequalities. So, 
in this respect our results do not give anything new. 

7.5. Our initial goal when starting the work on this paper was to 
construct potentials on Z'^, d > 3, such that N_{liav) = O^a^) but 
N^iHav) 7^ o(a2). In [18j this problem was mentioned as unsolved. 
We solved it here (in Section [H]) by using the sparse potentials. 

Sparse potentials were already used in the spectral theory of the 
multidimensional Schrodinger operator in the papers [I3l|llj. The au- 
thors studied the continuous spectrum of such operators, including the 
scattering for the pair (Ho,Hv). Their definition of sparseness always 
requires # {y G : \y\ < r} = o{r'^). Further requirements depend on 
the problem considered. The property (16.31) of almost orthogonality, 
which is the basis for our analysis, does not appear in their approach. 
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Note that the use of sparse potentials in the theory of one-dimensional 
operators has a long history, starting from the seminal work [TJ] by 
Pearson. 

7.6. The spectrum (including its discrete component) of the Hamil- 
tonian of the form (11. ip in dimension one is well studied, mostly due 
to the possibility to use the theory of Jacobi matrices. The first results 
on counting the number of eigenvalues outside the interval [0, 2] are 
due to Geronimo [6], [7], see also [9j, where a minor error in [7] was 
corrected. Other results in [H] concern the Lieb - Thirring type in- 
equalities, similar to those ioi d>3 described in Section [73] but much 
more elaborated. 

7.7. We do not touch upon the case d = 2 in this paper. Some facts 
can be derived from the corresponding results for the continuous case, 
by means of the technique developed in Section HI However, the results 
for d = 2 are incomplete, as well as in the continuous version of the 
theory, see the survey [18], and they leave many natural questions 
unanswered. In particular, the approach based upon sparse potentials 
does not work for = 2, so that in this case we still have no examples 
giving N{Ho,v) = 0{a) but N{Ho,v) 7^ o(a). 
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